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ABSTRACT 

We have used the visible integral-field replicable unit spectrograph prototype 
(VIRUS-P), a new integral field spectrograph, to study the spatially and spectrally 
resolved Lyman-a emission line structure in the radio galaxy B2 0902+34 at z = 3.4. 
We observe a halo of Lyman-a emission with a velocity dispersion of ^250 km 
extending to a radius of 50 kpc. A second feature is revealed in a spatially resolved 
region where the line profile shows blueshifted structure. This may be viewed as either 
HI absorption at ~-450 km s^^ or secondary emission at ~-900 km s~^ from the pri- 
mary peak. B2 0902+34 is also the only high redshift radio galaxy with a detection of 
21 cm absorption. Our new data, in combination with the 21 cm absorption, suggest 
two important and unexplained discrepancies. First, nowhere in the line profiles of the 
Lyman-a halo is the 21 cm absorber population evident. Second, the 21 cm absorption 
redshift is higher than the Lyman-a emission redshift. In an effort to explain these 
two traits, we have undertaken the first three dimensional Monte Carlo simulations of 
resonant scattering in radio galaxies. We have created a simple model with two pho- 
toionized cones embedded in a halo of neutral hydrogen. Lyman-a photons propagate 
from these cones through the optically thick HI halo until reaching the virial radius. 
Though simple, the model produces the features in the Lyman-a data and predicts the 
21 cm properties. To reach agreement between this model and the data, global infall 
of the HI is strictly necessary. The amount of gas necessary to match the model and 
data is surprisingly high, > IO^^Mq, an order of magnitude larger than the stellar mass. 
The collapsing structure and large gas mass lead us to interpret B2 0902+34 as a pro- 
togiant elliptical galaxy. This interpretation is a falsifiable alternative to the presence 
of extended HI shells ejected through feedback events such as starburst superwinds. An 
understanding of these gas features and a classification of this system's evolutionary 
state give unique observational evidence to the formation events in massive galaxies. 

Subject headings: galaxies: high-redshift — galaxies: individual (B2 0902+34) — galax- 
ies: formation — galaxies: evolution — galaxies: active — radiative transfer 
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Introduction 



Progress in understanding high redshift radio galaxies (HzRGs, where we arbitrarily consider 
high to mean z^2) will advance our picture of both massive galaxy formation and evolution. HzRG 
are thought to be the progenitors of massive local cluster giant elliptical galaxies. Their radio signal 
allows detection over a huge range of cosm ic time. We will only discuss a few revela i it tren ds of 
HzRGs as extensive summaries are given by iMcCarthyl (Il993l ) and lMiley &: De Breucld (|2008l ). Ra- 
dio galaxies are found t o be associated with so me of the largest stellar mass galaxies at all redshifts 



(jLilly &: Longaid Il984l : iPentericci et al.l l200ll ) . The environments of massive radio galaxies have 



also been progressively trace d to higher redshifts with radio galaxies typically residing in cluster 



environments at a l l redshifts (iLongair &: Seldneiill979l : iPrestage h Peacocklll988l : iHill k, Lillvlll991 



Carilli et al .1119971 : iBestI 120001 ). HzRGs often show extended Lyman-a emission up to and beyond 
their radio emission radii. The Lyman-a emission probes warm, ionized gas. Against this emission 
is often seen line profile structure which serves as a probe o f the neutral gas through scattering 



and absorptio n. The properties of Lyman-a blo bs fLABs) (IKeel et al 



1991 



1999 



Steidel et al 



Barkana &: Loebll2003l ) are 



good : 



Palunas et al.ll2004l ) and the gas halos of quasars (jHeckman et al 
morphologically similar to those in HzRGs. 

Cooling radiation may be an important source of extended Lyman-a radiation, and an indica- 



tion of ongoing galaxy fo rmation as infalling gas continues to build the galaxy's mass (jHaiman et al 



2OO0l : IPardal et al] boOlh . The most direct signature of a high redshift cooling flow is extended 
Lyman-a emission powered by the decrease in gravitational potential energy from infalling mate- 
rial. This is a challenging observation and interpre tation as othe r sourc es of Lyman-a generat i on are 
usually stronger than the cooling radiation, but iNilsson et al.l (j2006l ) and ISmith JarvisI (j2007l ) 
show promising cases. Low red shift cooling flows have a long history of in vestigation primarily 



Peterson et al. 



2001 



Tamura et al. 



through X-ray observations fe.g. iFabianI 1994 ). but space based observations (IKaastra et al 



2001 



2001 



) have given evidence against most cases of low redshift, 
strong cooling flows and diminished the importance of this mechanism for galaxy grow th. Most 
massive galaxy growth is now thought to occur through major mergers. iLi et al.l (j2007l ) show by 
a semi-analytic model, for example, th e likely growth history of a galaxy like the SDSS quasar 
J1148+5251 at z=6.42 (jFan et al.ll2003l ) where gas rich major mergers dominate as the source of 
the mass growth and a starburst phase predates a quasar phase. Observationally checking and 
challenging theories for the relative contributions of mergers and steady gas accretion for different 
galaxy types and epochs is an important goal. 

In terms of evolution, AGN feedback has been proposed ( Valageas h Silk 1999h as the mecha- 



nisin to produce the tight black hole /stellar mass correlation ( Magorrian et al.l 



2000 



1998 



— — II 7\ 

Ferrarese fc Mer ritt"20 0i3|) and the difference between the theoretical (jKaiserlll986l ) and the ob- 



Gebhardt et al 



served (jMarkevitchlll99S : Arnaud Evrardlll999l ) relation between cluster temperature and X-ray 
luminosity. Radio galaxies in particular and their powerful jets have been invoked as candidates for 
causing feedback (jinoue k. Sasakill200ll : iRawlings k. Jarvisll2004l : ICroton et al.ll2006l ) by depositing 
mechanical and thermal energy to their surroundings. A signature of outflowing gas at a velocity 
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exceeding the halo escape speed would be an important observation of high redshift feedback. 



B2 0902+34 (jLillvl Il988l ) is a powerful high redshift {z f» 3.39) radio galaxy of moderate 
projected radio diameter (40 kpc, 5") with a huge, diffuse Lyman-a halo lOOkpc) and is one of 
the most thoroughly studied HzRGs. Three observations make B2 0902+34 a unique window into 
HzRG halos as a class. First, 21 cm a bsorption from H I gas against the radio source continuuin 
has been observed a t Z2icm = 3.3962 (jUson et al.lll99ll : iBriggs et al.l 119931 : ICodv &: BraunI 120031 : 
Chandra et al.l |200J) , but any corresponding feature from this HI population has not shown in 
the Lyman-a emission line structure. This is the only HzRG for which 21 c m HI absorption 



has been repeatably ineasur ed with se ven oth e r HzR Gs giving null detections (jRottgering et al 



19991 : iDe Breuck et al.ll2003l ). Second, ICarillil (|l995l ) used radio spectral index and polarization 
mapping to constrain the approaching, northern radio jet inclination to the observer to between 
30° and 45°. This is a particularly low inclination angle for a radio galaxy, but it does offer very 
different lin es of sight to the two radio jets and any other structures that may follow the radio 
jets. Third, IVillar-Martm et al.l (l2007al ) have shown the emission line ratios in B2 0902+34 to be 
consistent with pure AGN photoionization through the Lyman-a, CIVA1549, and HeIIA1640 lines. 
A subsample of their HzRGs show line ratios that requires stellar photoionization. These systems 
would have peculiar and complicated source distributions for Lyman-a emission, but B2 0902+34 
appears to display the relatively well understood geometry of Lyman-a emission powered by an 



AGN. Until recently, longs 



lines ( Martin-Mirones et al 



often seen in other HzRGs ( 



it spe ctroscopy had detected only Lyman-a, CIV, and Hell emission 



19951) and shown n o sign of a bimodal Lyman-a line profile which is 



van Qjik et al.lll997l ). Deeper spectroscopy and a different slit positio n 



angle revealed a spatially resolved region with a bimodal Lyman-a line profile (jReuland et al.ll2007l ) , 
with most of the emission still from a simple and quiescent line profile. We will from here forward 
define a quiescent Lyman-a line as one with a <1 000 km s~^ t o indicate emission which does not 



display turbulent, radio jet dominated kinematics. iReuland et al.l (|2007l ) have also reported a weak 
[O II1A3727 detectio n in the galaxy's ce nter. Previous observations in deep Lyman-a narrowband 
(|Reuland et"al]l2003l ') and K broadband (|Eisenhardt Dickinsonlll992l l imaging and spectroscopy of 
the [OIII] 5007A line (|Eales et al.lll993l ) have led some authors to tentatively interpret B 2 0902+34 



as a protogalaxy based on its flat spectral energy distribution. iReuland et al.l (|2003l ) showed a 
diffuse Lyman-a morphology with two emission peaks roughly perpendicular to the radio axis. 
Unlike the two other HzRGs in their sample, B2 0902+34 did not display filamentary structure in 
Lyman-a. 

We present the first integral field spectroscopy (IFS) data of L yman-a emission in B2 0902+34. 

This is the fo urth study of HzRGs ha lo s using IFS after the w orks o flAdam et al.l (119971) . lyillar-Martm et al 
(|2007bl ). and iNesvadba et al.l (|2008l ). ILc Fevre et ahl (|l996l ) and Ivan Breukelen et al.l ^20od ) have 
used IFS to search for Lyman-a Emitters (LAEs) around HzRGs. With knowledge of the full spatial 
distribution of the bimodal Lyman-a line profile, we create a simulation for resonant scattering of 
the Lyman-a photons through a halo of HI. We find a model that reproduces most of the important 
spatial and spectral properties of our IFS observation and strictly requires a collapsing structure 
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with large HI mass. The natural explanatory power of this picture for other observations, such as 
the 21 cm absorption data, leads us to claim that B2 0902+34 is a protogalaxy still experiencing 
initial collapse. 

In ^we describe our observations. In ^we present our new resonant scattering model and 
analyze the impact of our model on 21 cm measurements. We show in §1] that our data is best 
fit through resonant scattering against an infalling halo. In we propose further discriminating 
observations and review. Throughout, we use a standard cosmology with Ho=71 km Mpc~^, 
Qm = 0.27, and JIa = 0.73. This gives look back times of 11.8 Gyr, angular scales of 7.5 kpc 
arcsec"^, and a luminosity distance of 30 Gpc for B2 0902+34 at z = 3.4. 



Observations 



On 2007 UT March 14 and 15 we obs erved B2 0902+3 4 with the visible integral-field replicable 
unit spectrograph prototype (VIRUS-P) (iHill et al.l 120081 ) on the 2.7m Harlan J. Smith Telescope 
at McDonald Observatory. VIRUS-P is a spectrograph with an integral field unit (IFU) of 247 
fibers in a hexagonal pattern close pack with a one third fill factor. Fed at f/3.65, the fibers have a 
diameter of 4.1" and a 107x107 □" field of view. Seeing was ~2", significantly smaller than the fiber 
diameter. The first night was nominally photometric, while the second was not. The instrument 
configuration gave us 236 working fibers spread over the field. Astrometry of the fiber positions has 
been calibrated to 1.0" rms from a fixed, offset guide field. Observations of three offset positions 
(called dithers) fill in the area. We took three exposures of 1200s at each of six subdither positions. 
A subdither as defined here is a second dithering pattern centered on the fiber intersections of the 
first dither set. We use this mode to decrease the astrometric uncertainty of spatially unresolved 
sources, mitigate local CCD cosmetics, and subsample extended emission. The wavelength range 
was 3400-5680A with a spectral resolution of 5.0A. Wavelength calibration was done with eleven 
arc lamp lines of Hg and Ne yielding a ax = 0.4A. Twilight flats were used to remove pixel-to-pixel 
variation, fiber-to-fiber variation, and the fiber illumination profile. We have used our own custom 
written software for all reductions. Notably, this reduction does not perform any resampling or 
linearization of the data. T he estimates of the sky spectrum were made from neighboring objectless 
fibers an d fit with bsp lines (jDierckxlll993l ) in a similar manner to optimal longslit sky subtraction 



methods ( Kelsonll2003l ). We estimate a spectrally unresolved 5a line flux limit of 7 x 10 erg s ^ 



cm ^ and a surface brightness limit of 24.0 AB mag arcsec ^ for wavelengths longer than 4000A 
when under photometric conditions. 

We detected Lyman-a emission from the low surface brightness halo of B2 0902+34 in flfteen 
different flber positions. We only include fibers that show emission above a signal-to-noise ratio of 
3, calculated from the noise and counts within the full width half maximum of a single Gaussian 
fit. We perform Monte Carlo line fitting where 100 realizations of our data sampled from normal 
distributions with our estimated errors produce distributions in our fit parameters from which are 
extracted central values and dispersions. All uncertainty quotes are 1 sigma estimates by this 
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method and similar to the uncertainties from the covariance matrices. All quotes of line width 
remove the instrument resolution in quadrature. Table [1] shows single Gaussian fits to each fiber's 
spectrum. Figure [1] shows the positions of the detections against a narrowband image we have 
ta ken and the radio dat a of ICarillil (|l995l ). We confirm the general line profile properties reported 



m 



Reuland et al.l (|2007l ) for the region covering their slit, but our observations span the entire 2D 



structure. The entire halo contains an emission peak at 5339.0 ± 2.0A with a FWHM of 600±90 
km which is well fit by a single Gaussian. This primary emission is mostly spatially symmetric 
and smoothly varying in brightness near the radio core emission although slightly extended toward 
the NE. At ~ 15% of the strength of the primary emission, a secondary emission peak can be fit 
as a Gaussian at 5324.5=bl.7A with a FWHM of 630=b270 km s~^. This secondary emission exists 
in three fibers on the SW side of the halo. Throughout this work, we will refer to any flux with 
a line center above 5330A as primary emission and any flux with a line center below 5330A as 
secondary emission. For the three fibers with S/N high enough for a more complex fit, we give fits 
in Table [2] where the fitting function consists of two independent positive Gaussians. Depending 
on interpretation, this feature may be due to discrete emission sites or scattering and absorption 
by neutral hydrogen. To aid comparison, we give fits in Table [3] with a single positive Gaussian 
and an absorbing Voigt function. Figure [3] shows the data and fits to select, representative spectra, 
and Figure S] shows the brightness maps of the two spectrally deco r npose d signatures as positive 



Gauss ians. Our data show that the longslit spectroscopy of iLillvl (|l988l ). iMartin-Mirones et al, 



(jl995l ) , and iReuland et al.l (|2007l ) are consistent though only the latter found this secondary, bluer 
emission because of the different slit position angles under their observations and the different flux 
limits. The two dimensional spatial constraints on the secondary emission afforded by our data 
will be crucial to a proper interpretation. We have searched for nearby LAEs in the surrounding 
fibers but have found none within the field. We do not make a statistically significant detection of 
continuum in B2 0902+34 which places a lower limit on the rest frame equivalent width from the 
six brightest fibers at > 67 A(5 a continuum limit) as shown in Figure [H We have not attempted 
a flux calibration due to the non-photometric conditions, so the aperture to aperture variation 
in surface brightness is overestimated by our data, and the flux limit between apertures is not 
constant. 



3. Interpretation 



A bimodal Lyman-a line proflle may be generated by several effects. First. ICooke et al.l (j2008l ) 
give imaging and spectroscopy of a Lyman Break Galaxy where the line profile is similar to that of 
B2 0902+34, which they interpret as being due to multiple merging galaxies even though the line 
widths are surprisingly narrow and the morphology undisturbed. This possibility in B2 0902+34 
cannot be e asily discarded . The secondary Lyman-a peak and a peak of continuum emission in 
Figure 8 of iReuland et al.l (j2003l ) appear somewhat coincident. However, our detection of the 
secondary emission peak spans three fibers, which covers a much larger area than the continuum 
signal. We speculate that this continuum may be stimulated or entrained material in the larger 
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Table 1. Individual fiber Lyman-a detections in B2 0902+34 



Unique position 


RA 


Dec 


Central A 


Dispersion 


Normalized 


S/N 


9 


Degrees of 


number 


(J2000) 


(J2000) 


(A) 


(A) 


Flux 






Freedom 


1 


9:5:30.35 


34:8:01.9 


5339.24±7.1 


11.5±9.7 




3.0 


10.8 


11 


2 


9:5:29.71 


34:8:01.8 


5338.96±0.9 


1.1±1.5 


0.19±0.07 


5.2 


6.4 


3 


3 


9:5:30.05 


34:7:55.1 


5337.24±1.0 


9.8±1.1 


l.OOitO.ll 


18.2 


25.0 


17 


4 


9:5:30.03 


34:7:59.7 


5338.60±0.6 


3.5±0.6 


0.39±0.05 


9.4 


1.1 


5 


5 


9:5:30.36 


34:7:53.0 


5338.84±2.2 


4.4±3.6 


0.40±0.20 


3.3 


2.5 


4 


6 


9:5:29.73 


34:7:53.0 


5332.51±2.8 


8.8±3.2 


0.33±0.11 


5.3 


13.1 


11 


7 


9:5:30.35 


34:7:57.6 


5338.45±1.0 


7.9±2.0 


0.68±0.13 


10.9 


11.3 


11 


8 


9:5:29.71 


34:7:57.5 


5339.27±1.5 


4.3±2.0 


0.26±0.07 


3.9 


7.7 


5 


9 


9:5:30.04 


34:7:50.8 


5337.66±1.3 


2.3±2.2 


0.15±0.11 


3.1 


13.9 


5 


10 


9:5:29.81 


34:7:59.7 


5337.92±1.0 


2.5±1.0 


0.20±0.05 


5.1 


4.4 


4 


11 


9:5:30.15 


34:7:53.0 


5337.66±1.3 


2.3±2.2 


0.25±0.09 


6.5 


5.7 


3 


12 


9:5:30.13 


34:8:01.9 


5336.14±5.6 


6.4±5.1 


0.22±0.13 


4.6 


16.0 


7 


13 


9:5:30.46 


34:7:55.2 


5339.12±2.2 


4.0±1.6 


0.19±0.07 


4.3 


3.8 


3 


14 


9:5:29.83 


34:7:55.1 


5338.03±2.9 


14.1±1.9 


0.68±0.16 


8.6 


39.8 


21 


15 


9:5:30.13 


34:7:57.6 


5338.96±1.0 


4.4±1.9 


0.33±0.09 


7.9 


1.7 


4 



Table 2. Select regions in B2 0902+34 fit with two Gaussians 



Unique position 
number 


Central Ai 

(A) 


Dispersioni 

(A) 


Normalized 
Fluxi 


Central A2 

(A) 


Dispersion2 Normalized 
(A) Flux2 


2 
5 
12 


5339.22±1.5 
5338.98±4.0 
5344.37±1.5 


6.5±2.5 
4.3±3.3 
4.7±3.0 


0.78d=0.27 
0.15±0.10 
0.31±0.13 


5323..34±7.7 
5326.97±3.2 
5325.49±1.0 


3+5.0 „ iq+0.28 
■'■■^-3.3 "■-'-^-0.19 

2.8±2.7 0.14±0.09 
2.9±1.0 0.22±0.06 


Table 3. 


Select regions in B2 0902+34 fit with Gaussian and Voigt profile 


Unique position 
number 


Central Ae 

(A) 


Dispersion 

(A) 


Normalized 
Flux 


Central Xa 

(A) 


b Column density 
(km/s) (l/cm^) 


2 
12 


5335.79±2.0 
5335.77±2.7 


8.9±1.6 
8.1±3.4 


1.11±0.39 
1.25±1.00 


5329.86±1.7 
5335.40±1.4 


liO_iio »-3XiU _8.3xiol3 
9,.t.+262 n 4yi(,14+3.8xl0l* 
^''='-266 ^-^XIU _2.4xlOl4 
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radio jet as invoked as a leading candidate to explain the alignment effect ([Begelman &: CiofB 



19891 : IChambers &: MilevI 119901 : lMcCarthvlll993l ) for radio galaxies. In the case of our supposition, 



the Lyman-a emission in this region originates around the rear radio jet and has a kinematic 
signature unaffected by peculiar motion of the continuum component. The large HI halo and 
scattering process we postulate would also cause a frequency redistribution which would wipe out 
any signature of the original kinematics, be they from a merging component or a rear photoionized 
cone. Also, the primary Lyman-a peak does not have a counterpart continuum peak, making it 
at least unlikely to be affiliated with a merger event. The second possibility is scattering from 
discrete regions of neutral hydrogen such as an outbound shell due to a starburst su perwind, and 



has been modeled and observed for LABs ( Ohvama et al 



2003: 



Wilman et al.ll2005l ) and Lyman 



Break Galaxies (jShaplev et al.ll2003l ). Third, a robust feature of optically thick resonant scattering 
in sta tic hydrogen clouds is a bimodal line profile straddling the rest wavelength (jUrbaniak &: Wolfe 
) . The addition of bulk, distributed velocity fields can then process the escaping spectra and be 



1981 



used to ma tch a wide range of observed, high redshift galaxies (jAhn &: Led 120021 : 1 Verhamme et al 



20071 . 12OO8I ). In these models, the source of Lyman-a emission is either a central point or a spherical 



distribution surrounded by a scattering medium. None of these pictures can seemingly explain the 
differing line profile seen between B2 0902-|-34's Lyman-a and 21cm line profiles. We submit 
a fourth mechanism, related to the third mechanism, to generate Lyman-a bimodality which is 
particularly suited to radio galaxies: two spatially segregated emission cones within a resonantly 
scattering halo. Resonant scattering effects have been shown through Lyman-a and Ha imaging to 
be important in low redshift objects (jHaves et al.ll2007l) and specifically with bulk velocity structure 
as the primary factor determining Lyman-a escape (|Kunth et al.lll998l ). 



3.1. Physical Assumptions and Constraints 



The basic geome try we investigate is shown in Figure [2j This scenario is based on a proposal in 
van Qjik et al.l (jl997l ) , although that work does not investigate the key effects of resonant scattering. 
The Lyman-a emission in HzRGs is often strongest along the radio axis i n a biconal geometry. This 
geometry has been advocated on AGN io nizing photori counting ground s (jMcCarthy van Breugel 



19891 ) and for AGN unification models ([Barthel 



198S 



Antonuccilll993l ). Surrounding the cones is 



a large halo of neutral hydrogen gas. Our model postulates an opening angl e of 90° for each cone . 
Matching the spatial location of the bimodal emission in our data and that of lReuland et al.l (|2007l ) , 
we model the projection of the emission cones axis on the sky at 60° east of north. It is not clear from 
the radio data exactly where the radio axis projects, but given the large deprojection uncertainty 
in such a low inclination system and the complex bending of the radio image, this geometry is 
reasonable. The ionization cone radius could be constrained by the Lyman-a luminosity, but 
since the filling factor is highly uncertain we keep this radius as an open parameter. We model 
the den sity profile as baryo ns at 17% of the total mass following Navarro- Frenk- White (NFW) 
profiles (INavarro et al.lll997l ) with a thermal core and concentration parameters of c = 5 following 
Dijkstra et al.l (|2006al ). hereafter DHS. The input velocity field is parameterized as a simple power 
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law with Vbuik = Vamp X Vmr X (^/'^iiir)'^ where a is the index, r is radius, r^r is the virial radius, 
Vyir is the virial velocity, and Vamp is a scaling factor of order one. The input parameters and 
acceptable ranges we explore in common to DHS are total halo mass, Vamp, and a. The halo is 
isothermal at lO^K, and the neutral fraction in the ionization cones is 10~^, both appropriate to 
the common ranges of density in photoionized regions. We model a completely smooth distribution 
of gas as a simplify ing assumption. In a very similar type of problem, low redshift cooling flows, 



Fabian et al.l ()1984l ) have shown how smooth flow assumptions adequately model infalling gas. The 



linear relation between the average number of photon scatterings before escape and the line center 
optical depth as shown in Appendix A ensures that gas clumping is irrelevant to our results. 



3.2. Resonant Scattering Model 



Given our new data of spatially resolved two dimensional Lyman-a emission, we have sufficient 
information to create a simple model for this system based on propagation of Lyman-a photons 
through an optic ally thick halo of HI with geometry as in Figure [2l Lyrnan-g radiative transfer 



is a mature field 


( Harrineton 


1973: 


Urbaniak & Wolfe 


1981: ] 


Meufeld 


1991; 


Loeb & Rvbickil 


1999; 


Ahn et al. 


2nnn. 


200] 




2002I: 


Zheng; & Miralda-Escude 12002: 


Richline 


200.4 


Hansen & Oh 


2006; 


Dijkstra et al. 


2006a 1 


2; 


Verhamme et al. 


200fil: Tasitsiomi 


200f 


: Laursen et al. 


2008 


) with both an- 



alytic and computational results for simple geometries. Broadly, the model will work by emitting 
Lyman-a photons in two photoionized cones roughly aligned with the radio axes. A red and blue 
emission line relative to systemic will emerge for each cone absent a velocity field in the scattering 
gas giving four peaks in the line profile. The far cone's emission will be weaker and further from 
the HzRGs systemic redshift. Global dynamics of infall will enhance the two blue emission peaks 
and suppress the two red emission peaks. Outflow will produce a mirror image spectrum about 
the systemic redshift and suppress the two blue emission peaks and enhance the two red emission 
peaks. The relative wavelengths between the stronger, narrower emission from the near cone and 
the weaker, broader emission from the rear cone can then discriminate between infall and outflow 
even if the systemic redshift is unknown. This method provides constraints on the HI mass and 
velocity field by measuring the bimodal line properties and the systemic redshift if available. 

The study of Lyman-a resonant scattering can be profitably attacked through Monte Carlo 
methods. We base our computations primarily on DHS, where DHS used Monte Carlo methods to 
attempt discrimination between energy sources and conditions in the data on LABs. Our version 
of the Monte Carlo resonant scattering code differs in three ways from their description. First, we 
do not include deuterium in our simulations. Although they show convincingly that deuterium is 
important for certain wavelengths and conditions, our very high optical depth and temperature 
environments make its contribution negligible and a needless drain on computing performance. 
Second, all of DHS's geometries are restricted to spherical symmetry. As we will invoke the two 
cones of the AGN as coaxial with the cones of Lyman-a emission, our models are without spherical 
symmetry. Our simulations therefore are a function of inclination angle and require the same 
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number of simulation photons in each inchnation angle bin as in an entire DHS simulation to 
reach an equal S/N. Third, DHS have a geometry that yields an analytic expression between 
optical depth and traversed physical distance for thin shells of constant density and velocity gas. 
The general density and velocity field conditions in our model mean that we must integrate this 
distance numerically. We use a root finding algorithm with a relative accuracy of 10~® on a fifth 
order Runge-Kutta code to calculate optical depths and sites of scatter. It is this generalization 
that enables us to pursue arbitrary density, temperature, and velocity fields. It will be possible 
to leverage this generalization by coupling our code to the more realistic conditions output by 
hydrodynamic codes in future work, but for now we use relatively ad hoc parameterizations. We 
used, in common with DHS, an acceleration scheme with Xcrit = 3 as described in their §3.1.1. Our 
simulation work is an extension and application of the powerful method of DHS to a different class 
of object, HzRGs. We have run our code for the four test cases shown in DHS's Figure 1 and found 
agreement as demonstrated in Appendix A. 

Our main differences in operation from the DHS models are the changes in HI density in the 
presence of ionization cones and AGN photoionization as the dominant Lyman-a energy source 
determining the emissivity profile. Lyman-a emission powered by a cooling flow and its associated 
uncertainties are unnecessary in B2 0902+34 with the strong evidence for AGN photoionization. 
Within the ionization cones, Lyman-a is generated randomly in position following the recombination 
rate. 

Our simulation ends when a photon leaves the virial radius. Further ra diative transfer effects 



in the IGM ar e beyond this work's scope, although they can be important ( Dijkstra et al.ll200 6b: 



Barkanall2004l ). We briefly discuss the qualitative impact that IGM pr ocessing by infalling n eutral 



hydrogen outside the virial radius could have to our model. Figure 2 of lDijkstra et al.l (|2006bl ) gives 
transmission curves at 2; = 3 for a variety of impact parameters and for cases with and without a 
nearby bright, ionizing quasar. Multiplying this transmission into a single peaked emission profile 
can produce a bimodal line profile, but this is not a likely explanation by itself for B2 0902+34. 
Only emission lines that are well centered around the systemic redshift can be modified by the IGM 
into a bimodal profile as the transmission minimum occurs slightly redward of systemic. All our 
models produce Lyman-a emission that is blueshifted enough that the IGM processing could not 
create a bimodal spectrum. The higher 21 cm redshifts also imply that the Lyman-a emission is not 
near the systemic redshift. We also do not observe the spatial dependence with impact parameter 
that IGM processing would create. For these reasons and scope, we neglect IGM processing. 

We must review a few of the broad findings in DHS to understand our models. Radiative 
escape from a static gas cloud will produce a bimodal spectrum. A photon has a much longer 
mean free path when it is at an off-resonance wavelength, which the thermal kicks in a scatter can 
cause. This situation will produce peaks redward and blueward of systemic with higher optical 
depths pushing the peaks out to more separated wavelengths. If the same situation has global 
infall (outflow) in the scattering gas, the photon moving outward and being scattered will have on 
average a blue-(red-)shift in the frame of the atom. A photon formerly in the red (blue) wing will 
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move into resonance and have its travelling distance supressed, while a photon formerly in the blue 
(red) wing will be pushed further into the scattering wing and travel further before its next scatter. 
In practice, the suppressed emission peak is likely to be far below detection limits for the range of 
optical depths we investigate. Therefore we return to a solitary emission line with no information 
remaining on the system's systemic redshift for cases of spherical symmetry. As DHS acknowledge, 
there is a degeneracy in observation between systems in infall and outflow. We turn to a stark 
physical difference between the DHS modeled LABs and our HzRGs to solve this problem. HzRGs 
frequently have two discrete regions of emission as ionization cones powered by the AGN. In a 
low inclination angle system such as B2 0902+34, the optical depth difference between these two 
regions to an observer can be substantial. Such a system should exhibit a weaker, spatially restricted 
emission profile caused by the rear ionization cone superimposed on a stronger, nearly circularly 
distributed emission profile caused by the near ionization cone. The relative order in wavelength 
between the weak and strong emission peaks will point in the direction of the unobserved, systemic 
redshift. So, HzRGs may be unique as high redshift halos where observation can distinguish without 
ambiguity between infall and outflow, solely from Lyman-a emission profiles. B2 0902+34 provides 
a further test of this idea since it is the only HzRG with a detection of 21 cm HI absorption. 



3.3. Evaluation of Model 

In Figure [3] we give the observed and simulated spectra in identical apertures for a particular 
model. There are five tunable parameters to our model: the virial radius, the ionization radius, the 
velocity strength, the power law index to the velocity field, and the photoionized cones' opening 
angle. We have coarsely run models over the likely halo parameter space, but we only present one 
model which captures many of the necessary observed traits. This model realization uses 700,000 
photons generated in a 6 x IO^^Mq halo which corresponds to a virial radius of 134 kpc and a 
virial velocity of 437 km s~^. The ionization radius extends to 101 kpc for this realization. We 
use an infall model where the velocity varies as the radius to the power a where here a = 0.5 and 
with a maximum radially inward velocity at the virial radius of 656 km s~^. The cone opening 
angle is 90°. We have one further parameter which we hold constant as it is well constrained by a 
previously discussed observation: the system's inclination angle. We believe 9i = 45° as most likely 
given the absence of an observed broad line region, so in practice we bin the spatial projections and 
wavelengths for all photons that exit within 35° to 55° of an ionization cone axis. For our model, 
we have assumed a systemic redshift of Zgystemic = 3.3990. We have done this after judging the 
velocity offset that the 21 cm absorption measurement should have in ^3.4[ The various redshifts 
are reviewed in Table HI 

This model reproduces many features of the observation. Most importantly, the size of the split 
in wavelength between the primary and secondary peak is recovered in the bottom panel of Figure 
[3l The relative intensities between the primary (near cone) and secondary (rear cone) emission 
are also fairly well matched. The regions of the halo not in projection of the receding radio cone 
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also show the proper behavior as demonstrated in the top plot of Figure [3j The emission there is 
not bimodal, and the peak is nearly constant at the same wavelength as the redder peak in the 
bimodal region. Several features in the model clearly do not meet the data, and our model is not 
to be considered the optimum match in parameter space but instead an early attempt. The model 
data all display a blueward skew, especially in Fiber 4 (top panel. Figure [3|) and regions away from 
either photoionized cone. Our data's S/N does not allow a useful measure of skew, but we predict 
that a much deeper and perhaps higher spectral resolution observation would show skew. This lack 
of skew is currently the most outstanding mismatch between our model and data. 

We consider our results as surface brightness distributions in Figure [H The left panel shows 
the observed decomposition between the primary and secondary peaks. The middle panel shows 
a similar decomposition for our model, where emission is decomposed spectrally. The right panel 
shows the model with the emission decomposed by the photoionized cone of origin. When we 
fairly decompose the model as in the middle panel, the region of secondary emission is pushed out 
slightly from the HzRG center and the relative intensity from secondary emission is underpredicted. 
Other models we have run have matched the spatial contours more closely, but have other failings 
such as underpredicting the bimodal wavelength separation, giving much broader lines, or showing 
even more skew. The modest S/N's so far observed for this dim object will inevitably cause 
confusion between skew, line width, and a bimodal line profile, so our unoptimized model parameter 
choices are appropriate to the state of certainty in the current data. The spatial distributions and 
relative intensities are fairly well matched at the level of the coarse spatial sampling allowed by our 
observation's fiber size. The offset spatial peak from the near ionized cone is also a robust match 
and explains why the Lyman-a peak and the radio core do not coincide. 

Other parameter choices have some important trends. We cannot use a much smaller total 
halo mass, because the rest frame wavelength separation in the bimodal regions drops from 4A to 
lA for a total halo mass of 3 X Larger halo masses are allowed and give more freedom 

in the other model parameters to reach the crucial rest frame wavelength separation of 4A, so we 
consider our gas mass to be a lower limit. Larger a values produce narrower emission lines with 
more skew, while a value of -0.5 creates lines of huge (40A) dispersion and spectral shifts (thousands 
of km s^^ from systemic) that are incompatible with the data. Smaller Ri to Ry ratios seem to also 
lower the rest frame 4A line profile separation unacceptably. A smaller value of 9o = 60° gives the 
appealing feature of more cleanly separating the near and rear photoionized cones spectrally much 
as the right panel of Figure H] shows, but the emission not in projection of the near cone begins 
to show line centers improperly centered around the wavelength of secondary emission contrary to 
observation. The most robust result of these models is that a rest frame line separation of 4A can be 
accomplished with resonant scattering from HzRG cones with a very large HI mass (lO^^M©), but 
a separation of much more (say, SA) can not without masses an order of magnitude larger. 
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3.4. 21 cm Predictions 

Any correct physical model of B2 0902+34 must take account of the 21 cm absorption data. 
Table [3] shows that fits of absorption to t he Lyman-a liii e prof ile do not match the absorbing 



population of HI evident in the 21 cm data. ICodv &: BraunI (|2003l ) detect an absorption feature of 



FWHM 120 km s"^ at z = 3.3962 with a column density of 3 x lO^^ cm under an assumed spin 
temperature of Ts=10^K. The relative velocities between the Lyman-a emission and the 21 cm 
absorbing HI and the absence of a signature from the 21 cm absorping gas in the Lyman-a profile 
both give strong constrains on possible geometries and conditions. 

Our model naturally reproduces the 21 cm absorption measurements without the need for a 
counterpart dip in the Lyman-a profile. By embedding the emission sites of Lyman-a in a larger 
neutral halo, the effects of scattering will not leave an absorption-like imprint on the Lyman-a 
profile, but they will do so on the 21 cm profile against the radio continuum emission. All existing 
21 cm data on B2 0902+34 are spatially unresolved. Figure [6] predicts the 21 cm maps that could 
be made using radio interferometry. The exact source distribution at 21 cm (323 MHz observe d 



frame) is unknown, but we estimate it to be the same as the 1.65 GHz distribution of lCarillil (|l995l ). 
We assume the radio source plane to be the sky plane rotated onto the radio axis at an inclination 
of 45°. We use the same model parameters as introduced in ^3.31 We find a column density of 
1.0x10^^ to l.Ox lO^^cm"^, a dispersion of 135 to 210 km s~^, and a velocity shift from systemic of 
-190 to -275 km s^^. These column densities would make this object a Damped Lyman-a system 
with a radius > 50 kpc. The highest intensity radio emission resides is in a region of large gradients 
for the absorption properties, so slight changes in source distribution or rotation in the plane of 
the sky cause large changes in the integrated values. We give the ranges that misalignments of 
±1" can create. The predicted line width is a factor of 2 above the measurement, but the rather 
unconstrained a parameter in the velocity field could likely remedy the mismatch as it takes a value 
closer to zero. Otherwise, the observations are bounded by these ranges, and the velocity shift tells 
us that our model self-consistently has both Lyman-a emission peaks blueward of both the 21 cm 
redshift and the systemic redshift. It is the velocity shift of 190 km for our fiducial astrometry 
that led us to the earlier assertion of Zgystemic = 3.3990. When building our resonant scattering 
model, we only let the 21 cm measurements guide our choice of systemic redshift. It is encouraging 
that our model designed to match the bimodal Lyman-a line profiles recovers so well the spatially 
unresolved 21 cm absorption signatures. 

Our most robust prediction comes from the southern radio lobe. Approximately 10% of the 
Table 4. Redshift estimates in B2 0902+34 

Primary Lyman-a Secondary Lyman-a 21 cm absorption Systemic from model 



3.3918 



3.3795 



3.3962 



3.3990 
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radio signal originates from this region in a ~l" diameter aperture judging from the 1.65GHz 
data. Assuming the absorption strength is uniform over the halo, this aperture should have an 
absorption strength of 1.6mJy. Spatially resolved observations would need an unattainably long 
1100 hours of exposure with 4 dishes from the Very Long Baseline Array (VLBA) and the phased 
Very Large A rray (VLA) to reach the same S/N of 6 as in the spatially unresolved absorption 
measurement ( Cody Sz BraunlboOS ) at a spectral resolution of 62.5kHz (restframe 58km s~^) around 
the 323.053MHz observed frequency according to the European VLBI Network calculatoi0. Adding 
more VLBA dishes or other facilities, such as Greenbank or Arecibo, does not lower the exposure 
time as the increased baselines overresolve the southern lobe. However, the brighter northern jet 
could be measured at a S/N of 5 in 40 hours under the same configuration, and we are pursuing that 
observation. A spatially resolved 21cm measurement will give important confirming or falsifying 
evidence for a vast and massive inf ailing gas halo. 



4. Implications: Infall Versus Outflow 



We are left with two serious, competing pictures: outflow from a 100 kpc HH region or infall 
of a 100 kpc HI region with embedded cones of HII. Deciding between these mechanisms has 
impact beyon d B2 0902+34. The b imodal line profile we see in B2 0902+34 is fairly common in 
HzRG halos. IVan Ojik et al.l (119971 ) found 11 of 18 such systems to carry the feature in longslit 
spectroscopy, although they all have not been studied for the confirming 2D spatial profiles as can 
be done with IFS. 

The outflow picture requires shells or clumps of HI at a large radius from the emission halo. 
This gas scatters a fraction of the line emission, and the blueshift in the absorption trough argues 
f or outflow. The en ergy mechanism might be a starburst-driven super wind as discus s ed ab ove. 



Wilman et al.l (|2005l ) has applied such an explanation to LAB2 as has iReuland et al.l (120071 ) for 
several HzRGs including B2 0902+34. If correct, these systems may be showing feedback in action 
with gas being heated and driven out. We see three reasons that this model is likely inappropriate 
to B2 0902+34. First, the 21 cm absorbing HI has no natural home in this picture. The radio 
emission, which extends to ~40 kpc must lie behind this bit of HI, but the Lyman-a photons must 
be generated in front or must at least have many scatterings aft er the absorptiq r i as no spectra shows 
Lyman-a absorption at the 21 cm redshift and column density. Reuland et al. (j2007l ) recognize this 
and predict a small, dense cloud of HI around the AGN. This too is unli kely, as the a mount of 
spectrally flat radio core emission available fo r absorption in t he da ta of ICarillil (jl995l ) and the 
amount required by the 21cm measurement of ICody &: BraunI ([2003) are modestly incompatible. 
Second, the surface brightness profile has the wrong slope for unscattered, extended emission. As 
shown in Figure [71 the observed surface brightness profile follows a power law. Taking an NFW 
profile for density and integrating the square of the density through the halo, assuming an optically 



http://www.evlbi.org/cgi-bin/EVNcalc 
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thin case and highly extended photoionization, we find surface brightness profiles that are too 
steep for all reasonable halo masses. Our simulation of resonant scattering gives a nearly power 
law profile with an appropriate slope. Third, the geometry of an absorbing cloud would have to 



be fin ely tuned to produce the observation. Based on our spectroscopy and that of iReuland et al 



( 20071 ). the primary peak's center stays very constant at 5339.0 it 2.0A in all positions. In the NE, 
there is no secondary emission peak so the absorption must be total and go quite blue. In the SW, 
the absorption wou l d hav e to be weaker, narrower, and redder than in the NE to produce Figure 
4 of iReuland et al.l ([20071). The outflow picture does not produce a natural explanation for the 
spatial variation in the line profile. 



Our result is surprising compared to the study of lNesvadba et al.l (|2008l ) with data from three 
other HzRGs. The authors find extended emission in rest-frame optical lines and kinematics in- 
dicative of outflow. A possibly crucial difference between the datasets is that for B2 0902-1-34 the 
emission line regions extend far beyond the radio lobes, while the other three HzRGs are bounded 
by radio emission and significantly elongated along the radio axis. If we speculate that both our 
and their conclusions are true, the next questions becomes whether the same HzRGs show different 
dynamical signatures through different lines or whether different HzRGs show different dynamical 
signatures across all their lines. 

Th e infall picture has been suggested before for HzRGs (e.g. [Humphrey et al I2OO7I ). IVillar-Martin et al 



(|2007bl ) produced a model of optically thin biconal emission in Lyman-a to recover their observa- 
tions of a peak emission position and peak velocity offset position. Their models, however, required 
postulation of a dense core that blocks emission from a rear cone that would have upset their trends, 
so their model is inappropriate to a bimodal line profile. They also could not produce lines at the 
width of their observations. We agree with the infall picture but use an optically thick regime of 
ionization cones embedded in one large, continuous neutral halo to explain all the available data on 
B2 0902-1-34 as we discussed above. If correct, this method has the capability to measure distributed 
hydrogen mass which we predict to be of the order lO^^Mf?) . This is 5 orders of ma gnitude higher 
than the warm and cool hydrogen populations estimates by lDe Breuck et al.l (120031 ) for the HzRG 
J2330-I-3927 from optically thin assumptions and Ly man-g line profile absorptio n specifically, and 
similarly high for the typical mass budget of HzRGs ([Miley &: De Breuckl l2008l ) . The neutral gas 
mass is dom inant by a factor of >17 compared to the estimat e of stellar mass in B 2 0902-1-34 of 
6 X 1O^°M0 dSevmour et al.ll2007l ). We note that in the work of lSevmour et al.l ([2003), B2 0902-^34 
has an usually low stellar mass compared to other radio galaxies at a similar epoch. Combined 
with our work here, this low stellar mass may indicate HzRGs exist in a fairly heterogeneous range 
of galaxy evolution, with B2 0902-1-34 at an earlier stage of evolution than the average HzRG. The 
unusually early evolutionary state and large HI halo for B2 0902-1-34 may also be why only B2 
0902-1-34 so far has shown 21 cm absorption. Our model also puts the system in a radically differ- 
ent evolutionary state from that predicted by outflow. Instead, we have infalling gas that is still 
building up what will become the galaxy's luminous mass. Star and dust formation have not yet 
built up to a dominant level to affect the Lyman-a line. The central AGN has not yet blown out or 
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heated a significant fraction of its environment in an act of feedback. The model mass would place 
B2 0902+34 as one of the larger halos that should collapse at the observed redshift. This picture 
does not deny HzRGs the role of feedback agents, but does show that radio emission can be output 
in a protogalactic phase as well and that AGN buildup may predate significant star formation. 
Further study on B2 09 02+34 may be important in deciding the origin of supermassive black holes 
(jDjorgovski et al.ll2008l ). and their coevolution with stella r bulges. Two likely paths of early black 



hole formation are as end biproducts of the earliest stars (iMadau Rees 



2001 



Ricotti fc Ostriker 



20041 ) or as gas that avoided mole cular cooling and f ragmentation and collapsed into a black hole 
without associated star formation (Silk &: Reeslll998l: iLoeb &: Rasidll994l : lEisenstein &: Loeblll995l : 



Bromm &: Loeb 



2003 



Koushiappas et al.ll2004l ). In the latter phase of galaxy evolution may 



exist where the stellar bulge and black hole masses are not tightly related. Although we do not 
have an estimate of B2 0902+34's black hole mass, the large radio continuum luminosity suggests 
an already large black hole. B2 0902+34, although we do not observe any associated cluster LAEs, 
may yet grow into a rich cluster as the normal cluster members begin their ow n star formation, 
possibly triggered by B2 0902+34 itself as described in lRawlings &: JarvisI (|2004l ) . 



While we advocate a resonantly scattering, infalling model to explain the extent of emission and 
the otherwise conflicting Lyman-a and 21cm absorption data, our model does require a minimum of 
10^^ M0 in HI subject to the modeled constraint of a maximum infall velocity of Vhuik 



"amp 



XV,, 



with here Vamp < 1.5. H aying this much HI inside a galaxy's v i rial radius conflicts with the 



classic galaxy formation scenario (jBinnevlll977l : iRees &: Ostriker 



1977 



Silklll977l ). In that picture, 



a shock near the virial radius should ionize incoming gas and heat it to ~ IQ^K where it stays in 
quasi-hydrostatic equilibrium. Especially in the most massive halos, the low density gas is expected 
to have a cooling time longer than the galaxy age. This scenario has been adj usted in recent years 
with rapidly cooling accretion flows showing up in semi-analytic simulations (jCroton et al.l 120061 ) 
and cold accretion in smoothe d particle hydrodynamic s (SPH) simulations where a virial shock 
is mass conditionally unstable (jBirnboim &: Dekell |2003| ) and filaments provide many paths where 



cold gas (< 10^ -RT) can f low to the formiri g galaxy's center ( Fardal et al 



2001 



Ocvirk et al 



Keres et al. 2005 



200. 



Th e abil ity of 



Dekel k BirnboimI l2006l : iDekel et al.l l2008l : iKeres et al.l l2008l : 

cooling and cold flows to exist have a red shift and hal o mass dependence. ICroton et al.l (|2006l ) give 
Mhaio < 2 - 3 X lO^^M© for z=0-6 and IXeres et al.l (|2005l 'l give Mhaio < lO^^-^M© at z=0 with 
some redshift evolution. In either case, our advocacy for a Mhaio > 5 x lO^^M© at z=3.4 appears 
to be in conflict with theory. Instead, the growth mechanism of such a massive halo is expected to 
be mergers. Future radiative transfer simulations of Lyman-a in more realistic, filament condition s 



may alleviate this mismatch in mass scale, although since the velocity fields in Keres et al. (j2005l ) 
are of the same order as those we modeled, an easy solution is not obvious. The potential conflict 
with theory at such a large mass makes B2 0902+34 a compelling target for further observations 
and simulations. 
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5. Conclusions 

We have made new spatially resolved IFS observations of Lyman-a in B2 0902+34. We find a 
spatially resolved region of weaker and bluer emission superimposed on a symmetric halo of primary 
emission. We have shown that all properties of B2 0902+34 can best be fit with a resonant scattering 
infall model with a large, mostly neutral, hydrogen mass. The sign of the velocity field is relatively 
robust and model independent. Signal seen through higher optical depths becomes increasingly 
bluer indicative of infall and resonant scattering. We match our data with a model dependent HI 
mass of > IO^^Mq which is >17x larger than the stellar mass. For three reasons an outflowing 
HI scenario cannot fit our observations: 1) the strength of the observed 21cm absorption and the 
Lyman-a profile are incompatible with a single HI population if optically thin and the velocity 
offsets between the 21cm signal and the higher optical depth Lyman-a signal carry the wrong 
sign for outflow if optically thick, 2) the surface brightness profile, with the assumption that the 
gas follows the dark matter, would be steeper than observed if optically thin, and 3) the bimodal 
Lyman-a line profile exists only in a small, non-central region of the galaxy which is incompatible 
with an outflowing shell geometry. The 21 cm HI absorption has been a powerful confirmation of the 
kinematics through its redshift offset from the Lyman-a halo emission, and we believe observational 
progress in understanding this class of object will benefit from further deep 21 cm measurements on 
other HzRGs. This new scenario of a massive, infalling HI halo with resonant scattering places B2 
0902+34 in a much earlier evolutionary state than the previously formulated starburst superwind 
picture, and we classify B2 0902+34 as a protogiant elliptical galaxy in an early collapse phase 
before significant feedback has disrupted star formation and AGN fueling. Observations of other 
HzRGs with IFS will be important in correctly classifying systems between formation and feedback 
modes and allow detailed observational tests of these two processes. We are observing a large 
sample of HzRGs with VIRUS-P in order to undertake a similar analysis on other systems. 

Further types of data on this system will allow confirmation of our infalling, resona nt scattering 



model. Polarization measurements of Lyman-a could confirm the scattering process (jLee &: Ahn 



19981). The bright regions of HzRGs have been found to have modest levels of pol arization (^10% 



0, 



Vernet et al.l (|200ll )) while scattering simulations in similar situations to our model (jPiikstra Loeb 
20081 ) predict higher (~40%) values. A measurement between HzRGs with and without emission 
line halos may demonstrate this split. The key, potentially falsifying observation for this resonant 
scattering, infalling model of B2 0902+34 will be radio interfero metry spectra l imaging to spatially 



resolve the distribution of HI gas through its 21 cm absorption. ICarillil (119951 ) discussed the utility 
of such a measurement, and we have now given an observationally motivated model that predicts 
spatially extended HI. 



This material is based upon work supported under a National Science Foundation Graduate 
Research Fellowship. This work is supported in part by Texas Advanced Research Project Grant 
number 003658-0005-2006. The simulations used in this work were carried out at the Texas Ad- 
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Facilities: Smith (VIRUS-P). 



A. Tests of the resonant scattering code 



To provide compatibility between our work and previo us authors', we r un the same set of four 
tests on our 3D Monte Carlo resonant scattering code as iDiikstra et alj (|2006al ). The first test 
simulates the emergent spectra for centrally created Lyman-a photons in a uniform density spheri- 
cally symmetric static halo of HI with various line center edge-to-center optical depths at T=10K. 
We disabled the acceleration scheme for tests 1-3. The second test calculates the redistribution 
function for single scattering events, which is the probability distribution function for an output 
frequency given an input frequency and gas temperature. The third test finds the mean number 
of scatteri ngs prior to escape for a slab of gas of various line center optical depths. The weighting 
scheme of I A very Sz Housd () 19681 ) was used for test 3 where the number of accumulated scatterings 
and the probability of escape from each scattering site is retained as opposed to solely keeping the 
properties of each Monte Carlo run's escaping photon. The fourth test gives the e mergent spectrum 
from an infinitely large object under Hubble exp ansion as origi n ally sim ulated in lLoeb &: Rvbicki 
(jl999l ). More details and references are given in IDiikstra et al.l (|2006al ) Appendix B. We give the 
four tests in Figure [51 
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Fig. 1. — The fibers with detected Lyman-a emission are plotted over a narrowband Lyman-a 
image of the radio galax y cente r . For reference, log scale contours are drawn of the 1.65 GHz image 
at 0.15" resolution from lCarillil (|l995l ) with range 0.36-92.16 mJy/beam. Each fiber has a diameter 
of 4.2". 
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Fig. 2. — The geometry of our simulation. 6i is the inchnation as constrained by ICarillil (|l995l ). 

is the system's virial radius. Ri is the ionization radius of the cones. 9o is the opening angle of 
the ionization cones assumed here to be 90°. Ri, Ry, 9o, and two variables controlling the velocity 
field are the model's five tunable parameters. 
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Fig. 3. — The spectra of both our IFS data and resonant scattering simulation from select fibers. 
The thin dotted line indicates the redshift measured for 21 cm absorption. Solid Observed. Thick 
dotted Simulated. Top Fiber 4, representative of most of the halo area where neither cone directly 
projects. Middle Fiber 7, representative of regions spanned by the near photoionized cone. Bottom 
left Fiber 14, representative of regions spanned by the rear photoionized cone. Only these regions 
show a bimodal line profile. Bottom right A region in the simulation near fiber 14. While we did 
not have a fiber on this exact position, the bimodal spectrum here is a good match to its neighbor 
in fiber 14. 
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(a) Observed (b) Simulated, Spectral Decomposi- (c) Simulated, Origin Cone Decompo- 

tion sition 

Fig. 4. — Log contours of surface brightness. Left Fits to the data. The primary and secondary 
peaks are respectively dark and hght. Middle The simulations described in ^ where the same type 
of spectral decomposition as for the data between primary and secondary is performed. Right The 
same simulations, but now with a decomposition regardless of the spectral signature but where the 
dark contour follows emission from the near photoionized cone and the light contour from the rear 
photoionized cone. This optimistically represents a clean detangling of spectral signal that a model 
with more optimized parameters may accomplish. 
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Fig. 5. — The joint detection of Lyman-a emission in a subset of the central fibers. The equivalent 
width limit is shown. Each row in the figure's bottom half corresponds to the spectrum from a 
fiber with a unique position on sky. In the labeling of Figure [H these fibers are numbers 4, 15, 3, 
7, 14, and 11 from top to bottom. 



Fig. 6.— The predicted 21 cm maps for B2 0902+34. The 1.65 GHz data of lCarihil (119951 ) is plotted 
as the dashed contour as a hkely source distribution. Top Left Column density where the log of 
N(HI) in cm~^ is given. Top Right Velocity dispersion in km s^^. Bottom Velocity shift in km s~^. 
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Fig. 7. — Observed and simulated surface brightness profiles. All models have been scaled in 
intensity by least squares to the data. The optically thin integrals through NFW halos fail, though 
the resonant scattering model reproduces the observed, shallower profile. 
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Fig. 8. — The four tests of our 3D Monte Carfo resonant scatter code. Top left The emergent 
Lyman-a spectra from uniform density, static halos at T=10K with various Hue center edge-to- 
ce nter optical de p ths. T he histograms are from our simulations. The lines are from function 9 



m 



Dijkstra et al.l (|2006al ). Top right The redistribution function for Lyman-a scattering with a 
ga s temperature o f T=10fK. The histogram is our code output while the curves are from theory 
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Hummer! (jl962l ): iLed ()1974l ). Bottom left The number of scatterings necessary for a Lyman-a 



photon to es cape a neutral, s tatic slab of HL The points are from our simulations. The line is 
from theory (jHarringtonlll973l ). Bottom right The emergent Lyman-a spectrum from an infinite, 
uniform density HI halo under Hubble expansion at z=10 with T=10K. 



